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ABSTRACT 



The high-frequency-peaked BL-Lacertae object lES 0806+524, at redshift 
z=0.138, was observed in the very-high-energy (VHE) gamma-ray regime by 
VERITAS between November 2006 and April 2008. These data encompass the 
two-, and three-telescope commissioning phases, as well as observations with the 
full four-telescope array. lES 0806+524 is detected with a statistical significance 
of 6.3 standard deviations from 245 excess events. Little or no measurable vari- 
ability on monthly time scales is found. The photon spectrum for the period 
November 2007 to April 2008 can be characterized by a power law with photon 
index 3.6 ± l.Ostat ± 0.3sys between ~300 GeV and ~700 GeV. The integral flux 
above 300 GeV is (2.2 ± 0.5stat ± 0.4sys) x 10"^^ cm~^ s~^ which corresponds to 
1.8% of the Crab Nebula flux. Non-contemporaneous multiwavelength observa- 
tions are combined with the VHE data to produce a broadband spectral energy 
distribution that can be reasonably described using a synchrotron-self-Compton 
model. 

Subject headings: BL Lacertae objects: individual (lES 0806+524) — gamma 
rays: observations — X-rays: galaxies — ultraviolet: galaxies 



1. Introduction 



Active galactic nuclei (AGN) are galactic cores with a luminosity that dominates the host 
galaxy. These central engines are believed to be powered by superm assive black holes which 
draw surrounding matter into an accretion disc. The unified model (jUrry fc Padovanilll995l ) 
of AGN implies that the observed emission from AGN is a strong function of the observation 
angle in relation to the relativistic jet, with the highest energy observed emission beamed 
along the direction of the jet. 
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lES 0806+524 was identified as a BL Lacerta e object (ISchachter et al.lll993[) using radio 
observations from tlie Green Ba nk 91-m telescope (IBecker et al.lll99ll ) and X-ray observations 
from tfie Einstein Slew Survey (lElvis et al.lll992l ). The redshift of the host galaxy is z=0.138 
( iBade et al.lll998l ). The blazar lES 0806+524 was suggested as a candidate VHE gamma- 
ray source based on the presen ce of both high-energy electrons and sufficient seed photons 
( ICostamante fc Ghisellinill2002l ). That work predicts an intrinsic flux in the VHE regime of 



E>0.3 TeV 



1.36 X 10-11 cm-2 s-i. 



Very high- energy gamma-ray observations of lES 0806+5 24 reported by the Whipple 
Collaboration (IHoran et al.l l2004l : Ide la Calle Perez et al.l |2003| ) indicated flux upper limits 
of ^E>o.3TeV < 1-4 X cyot'^ -sT^ and J^E>o.3TeV < 16.8 X 10-11 cm-^s"! from the 1996 and 
2000 observing seasons and ^E>o.3TeV < 1-47 x lO^^cm-^s-i from the 2000 to 2002 observ- 



ing seasons. Observations reported by the HEGRA Collaboration (lAharonian et al.l l2004l ) 
indicated an upper limit of J?-e>i.o9 



TeV < 43 X 10 11 cm ^ s 1 based on hour of observations. 



Obser vations of lES 0806+524 in late 2006 and early 2007 with VERITAS fICogan et al. 



2007al ) yielded evidence for weak emission with a statis t ical sig nificance of 2.5 standard de- 
viations in 35 hours. The data reported in ICogan et al.l (l2007al ) are combined with the data 
from the 2007/2008 observing season here. 



2. Observations and analysis 



VERITAS ( IWeekes et al.ll2002l ) is an array of four 12-m diameter imaging Cherenkov 
telescopes located at the Fred Lawrence Whipple Observatory in southern Arizona (1268 m 
a.s.L, 31°40'30"N, 110°57'07"W). 

Observations of lES 0806+524 were made between November 2006 and February 2007 
during the construction phase of VERITAS and from November 2007 to April 2008 during 



the f ull four-telescope operation. All observations were taken in wobble mode (IFomin et al. 



19941 ). where the target is offset from the center of the field of view by ±0.3° or ±0.5°. 
After quality selection to remove data suffering from poor weather or technical problems, a 
total of 65 hours of observations are available. These data comprise 5 hours with a T1/T2 
configuration, 25 hours with a T1/T2/T3 configuration, and 35 hours with a T1/T2/T3/T4 
configuration. Th ese data were analyz ed using the standard offline VERITAS data analysis 
package VEGAS (ICogan et al.l l2007bl ). Consistent results are obtained using independent 
analysis packages. 

Shower images of extensive air showers in each camera are gain-cor rected and cleaned 
(ICogan et al.ll2007d ) before being parameterised using a moment analysis (lHillaslll985l ). Prior 
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to stereoscopic reconstruction, a pre-selection retaining images with size greater than 75 
photoelectrons, with greater than 4 pixels surviving cleaning, and with a distance from the 
center of the field of view to the image centroid of less than 1.43° is applied. Events where only 
images from the two closest-spaced telescopes survive the pre-selection are rejected. Each 
shower is parameterise d using mean-scaled width (MSW) and mean-scaled length (MSL) 
( iKonopelko et al.lll995l ). Gamma-ray selection criteria of MSW < 1.16 and MSL < 1.36 are 
used to reject most of the background while retaining a large portion of the signal. These 
selection criteria were optimised assuming a source strength of 3% of the Crab Nebula flux. 



Background estimation is performed using the refiected-region background model (iBerge et al. 



20071 ). In this scheme, an integration region is placed around lES 0806+524, with back- 
ground integration regions distributed around the field of view. The number of events in the 
integration region is termed ON, the number of events in the background region is termed 
OFF, and the ratio of the integration areas is termed a. As the data comprise obs ervations 
with d ifferent abso lute wobble offse t s, a g eneralized version of equation 17 from iLi fc Ma 
( 1l983l ). as noted in lAharonian et al.l (120041 ). is used to compute the statistical significance of 
the excess. 

The entire VERITAS dataset from 2006 through 2008 reveals the existence of very- 
high-energy gamma-ray emission from lES 0806+524 with a statistical significance of 6.3 
standard deviations above an energy threshold of ~300 GeV. The total number of ON counts 
is 1543, the total number of OFF counts is 13311 and the effective ratio of the integration 
regions (accounting for the different wobble offsets) is a = 0.0975, yielding an excess of 245 
events. 

A distribution of the squared angular distance (referred to as 6^^) between the position of 
lES 0806+524 and the reconstructed shower direction is shown in Figure [H The integration 
region used in the refiected-region analysis is indicated by the vertical line situated at (P' = 

2 

0.013° . The scaled background counts are shown for comparison, and are computed using the 
centers of the reflected regions as the reference points. The distribution of excess gamma-ray 
counts is consistent with that given by a point source located at right ascension 8^9™59^ ± 9*^ 
and declination 52°19' ± 2'. There is an added combined systematic angular uncertainty 
of ~100". This is consistent, within error, wi th the position of t he radio source at right 
ascension 8^9™49.2^ and declination 52°18'58" JKovalev et al.ll2007h . 



Only those data taken with four operating telescopes are used for a spectral analysis, 
owing to the superior energy resolution afforded. The subset of data, taken during the 
2007/2008 observing period, reveals a total of 17 6 ON events, 13 34 OFF events with statis- 
tical significance of 4.4 a using equation 17 from iLi fc Ma Jl983h . where a = 0.09, yielding 
an excess of 55 events. The differential photon spectrum can be fit using a power law of 
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the form dN/dE = Fq x (E/400GeV)"' where Fq = (6.8 ± 1.7stat ± l-S^ys) x lO-^^^m-^s-^ 
and r = 3.6 ± l.Ostat ± O.Sgys. The fit yields a x^/dof of 0.07/2, where dof is the number of 
degrees of freedom. Upper hmits for t wo spectral p oints above ~700 GeV are calculated at 
the 90% confidence level according to iHelend (Il984j ). 



Absorption on the infrared compone nt of the extragalactic b ackground light results in 
an attenuati on of high-energy p hoton s (iGould fc Schrederl 119671 ). The absorption model 
according to iFranceschini et al.l (120081 ) is used to calculate the deabsorbed spectrum. This 
is achieved by scahng each fiux point according to = ^ohs^'^^^''^\ where E is the energy 
at that flux point, z is the redshift and r is a function describing the absorption model. A 
power law fit to the resultant fiux points results in a spectral index of —2.8 ± O.Sstat- 

The integral light curve above 300 GeV is shown in Figure [3] with the data binned per 
month. A constant fit yields a x^/dof of 6.78/5 indicating little or no measurable variability. 



3. Analysis of Swift UVOT and XRT Data 



Following the announcement of the discovery of lES 0806+524 by VERITAS in the 
very-high energy ga mma-ray regime (ISwordyll2008al ). lES 0806-f524 was observed by Swift 
( IGehrels et al.ll2004l ) in the ultraviolet to optical and X-ray energy bands using the UVOT 
and XRT between 2008 March 8 and 2008 March 12. Data reduction is carried out with 
the HEAsoft 6.4 package. The Swift XRT data were taken in photon counting mode at a 
rate below 0.6 counts/s, so photon pileup is not evident. The XRTPIPELINE tool is used 
to calibrate and clean the XRT event files. Source and background counts are extracted 
from circular regions of radius 30 and 40 pixels, respectively. The 0.6 to 10 keV energy 
spectra are fit by a power law with fixed Galactic column density of N h = 4.4 ■ 10^° cm~^ 
( Dickey &: Lockmanlll990l ). Marginal fiux variability is seen from 2008 March 8 to 12, with 
the photon index hardening from F = 2.67 ± 0.08 to F = 2.53 ± 0.07. 

Swift UVOT data is reduced with the UVOTSOURCE tool to extract counts, correct 
for coincidence losses, apply background subtraction, and calculate the source fiux. The 
standard 5 arcsec radius source aperture is used, with a 20 arcsec ba ckground re g ion. T he 
source fiuxes are dereddened using the interstellar extinction curve in iFitzpatrickl (119991 ). 
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4. Discussion and Modeling 

Non-contemporaneous data from the Swift UVOT0 and XRT are combined with the 
VERITAS data and an archival optical data point (from the Tuorla Im telescop^ to pro- 
duce a broadband spectral energy distribution in a vFj^ representation as shown in Figure 
HJ The Tuorla data were taken with an R-band filter and host-galaxy subtraction is applied. 
The spectral energy distribution ( SEP) is modeled using the equilibrium version of the one- 



zone jet radiation transfer code of iBottcher fc Chiang] (120021 ). with parameters appropriate 



for a pure synchrotron-self- Compton (SSC) model. This model assumes a population of 
ultrarelativistic nonthermal electrons (and positrons) injected into a spherical comoving vol- 
ume (referred to as the hloh). The injection spectrum is characterized by an injection power 
Linj {t) with an unbroken power-law distribution between energies 71 and 72 with index q. 
The jet moves with relativistic speed /3c and it is charaterized by Doppler factor D which is 
dependent on the line of sight. The injected particles suffer radiative losses via synchrotron 
emission and Compton upscattering of synchrotron photons. 

The SSC model has been fitted to the data, with different parameters obtained for 
the 2008 March 8 and 2008 March 12 data sets. For 2008 March 8, we obtain Li^j (t) = 
1.9 X 10*^^ erg s^^ with an injection spectrum described by 71 = 1.77 x 10^, 72 = 2 x 10^, 
spectral index q = 3.1 and a magnetic field strength of 0.39 Gauss. The corresponding 
parameters for the 2008 March 12 dataset are 1.6 x 10'^^ erg s~\ 71 = 1.6 x 10"^, 72 = 2 x 10^, 
injection spectral index q = 2.7 and magnetic field strength of 0.5 Gauss. The blob radius 
is 5 X 10^^ cm in both cases. A Lorentz factor of F = 20 is used, with the simplifying 
assumption that the Doppler factor D is equal to the Lorentz factor. The synchrotron peak 
is located at ~8.3xl0^^ Hz and the inverse- Compton peak is located at ~3.5xl0^^ Hz. The 
system is within a factor ~2 of equipartition in both ca ses. Absorption on t he ex tragalactic 



background light using the EBL scenario described by iFranceschini et al.l (120081 ) has been 
accounted for in the model. 



5. Conclusions 

VERITAS has observed the blazar lES 0806+524 for a total 65 hours from November 
2006 to April 2008 resulting in the discovery of very-high-energy gamma rays with a statistical 
significance of 6.3 a. The differential energy spectrum between ~300 GeV and ~700 GeV 



^Note that the host galaxy contribution is not important in UV. 
^http:/ /users. utu.fi/kani/lm/ 
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3.6 ± l.Ostat ± 0.3 



The 

X 10~^^ cm'^ which corresponds to 



Jsys- 



can be fitted by a relatively soft power law with index T 
integral fiux above 300 GeV is (2.2 ± O.Sstat ± 0.4sys 
(1.8 ± 0.5)% of the Crab Nebula fiux as measured by VERITAS above 300 GeV. Assuming 
absorption on the infrar ed component of the extragalactic background light according to 



Franceschini et al.l (120081 ). the intrinsic integral fiux above 300 GeV is (4.4±0.6stat ±0.5 



10~^^cm~^s~^ which is approxim ately one order of magnitude less than the fiux predicted by 
Costamante fc Ghisellinil (120021 ). The broadband spectral energy distribution can be fitted 
using a one- zone SSC model with standard parameters. 

Observations by the EGRET gamma-ray space telescope reveal gamma-ray emission 
from the region around lES 0806+524, however the large error box indicates that the emis- 
sion could be associa ted with either B 0803+5126 (z=1.14) or lES 0806+524, with the 



former being favored (ISowards-Emmerd et al.ll2003l ). Observations of this region with the 
Large Area Telescope on the Fermi Gamma- Ray Space Telescope should resolve the source(s) 
of gamma-ray emission in this region. 

VERITAS is the most sensitive instrument of its kind in the Northern Hemisphere, and 
it is ideally suited to observations of extragalactic objects. This h as been demons tr ated by 
the detection of three new BL Lac objects by VERITAS in 2008 JSwordvl J2008ah . Is wordy! 
(l2008bh . lAcciari et all J2OO 
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Fig. 1. — Distribution of the parameter the squared angular distance between the 
reconstructed shower direction and the location of lES 0806+524. The vertical line, located 
at ^2 = 0.013° indicates the size of the integration region. The simulated line indicates the 
expected shape of the 6*^ distribution for a point source of this strength. 
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Fig. 2. — Differential photon spectrum of lES 0806+524. The spectrum is well fit by a power 
law with index 3.6 ± l.Ostat ± O.Sstat- The dea bsorbed spectrum i s calcu lated by applying the 
extragalactic absorption model according to iFranceschini et al.l (120081 ). 
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Fig. 3. — Month-by-month integral flux above 300 GeV from observations of lES 0806-1-524. 
The chi-square probabihty of the straight hne fit is 0.24. 
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Fig. 4. — Broadband spectral energy distribution of lES 0806+524. The Swift data are taken 
on two separate nights, with a lower flux on 2008 March 8 (gray points) and a higher flux on 
2008 March 12 (black points). The solid and dashed curves are SSC fits to the 2008 March 
8 and 2008 March 12 Swift data, respectively. 



